Constraints on Kinetically Modified Inflation from WMAP5 
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Single field inflationary models with a non-minimal kinetic term (also called k-inflationary mod- 
els) can be characterised by the so-called sound flow functions, which complete the usual Hubble 
flow hierarchy. These parameters appear in the primordial power spectra of cosmological pertur- 
bations at leading order and, therefore, affect the resulting Cosmic Microwave Background (CMB) 
anisotropies. Using the fifth year Wilkinson Microwave Anisotropy Probe (WMAP5) data, we derive 
the marginalised posterior probability distributions for both the sound and Hubble flow parameters. 
In contrast to the standard situation, these parameters remain separately unbounded, and notably 
there is no longer any upper limit on ei, the first Hubble flow function. Only special combinations 
of these parameters, corresponding to the spectral index and tensor-to-scalar ratio observables, are 
actually constrained by the data. The energy scale of k-inflation is nevertheless limited from above 
to //inf < 6 X 10~®m,pj at two-sigma level. Moreover, for the sub-class of Dirac-Born-Infeld models, 
by considering the non-gaussianity bounds on the sound speed, we find a weak limit ei < 0.08 at 
95% confidence level. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

Accelerated expansion of the early Universe in scalar 
field models is usually achieved by domination of the 
field's potential over its kinetic energy. However, it has 
been shown that inflation can also be supported through 
modifications of the kinetic term, these models being 
known as "k-infiation" scenarios [l[ . Assuming the grav- 
ity sector to be described by General Relativity [the sig- 
nature of the metric tensor being (— , +, +, +)], the action 
describing such a "k-infiaton" ip{x^) reads 



scenarios that emerges from string theory. This is the 
case of the so-called Dirac-Born-Infeld (DBI) scenarios 
of brane inflation, which have recently received particular 
attention as effective field theory realizations of inflation 
within string theory [1 i, i, i, S, 0, i, H • The DBI infla- 
ton typically represents the distance (i.e., an open string 
mode) between branes in a higher-dimensional bulk, and 
the Lagrangian functional is given by the determinant of 
the induced metric along the moving brane, 



S" = — / d*x^/^ [R + 2kP {X, if)] 



P°-\X,f) = -T{^)Jl- 



2X 



+ T{^)-V{^). (3) 



(1) 



where k = Stt/to^ 
ing expression 



and where X is defined by the foUow- 



X 



-g'"'df,ipd^ip. 



(2) 



The quantity P{X^ ip) is any acceptable functional (see 
below) for the scalar field Lagrangian. For instance, the 
canonical scalar field model is recovered for P{X, ip) = 
X — V{ip), where V{ip) is the scalar field potential. 

The theory of k-inflation is interesting not only because 
it enlarges the class of inflationary models that one can 
consider but also because this is precisely this type of 



Here, T(ip) is a warping function encoding the bulk ge- 
ometry, and V{ip) is the interaction potential between the 
branes, assembled from contributions of various origin. 

In general, a model of k-inflation is sensible if the 
function P{X, ip) satisfies the two following conditions. 
Firstly, Px > 0, coming from the requirement that 
the Hamiltonian is bounded from below and, secondly, 
Px + '^XPxx > 0, which is necessary in order for 
the Cauchy problem to be well-posed, see for instance 
Ref. [1^] . Let us notice that the DBI models satisfy these 
requirements. Moreover, these two conditions play an 
important role at the perturbative level. Indeed, all ki- 
netically modified inflationary models share the common 
feature that the "speed of sound" 
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(4) 
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is generically different from the speed of light and may be 
space-time dependent. Clearly, the two above-mentioned 



inequalities guarantee that c^ > 0. This is relevant for 
the treatment of cosmological perturbations. Indeed, for 
a flat Friedmann-Lemaitre-Robertson-Walker (FLRW) 
universe, it was shown in Ref. [ll|, that the Mukhanov- 
Sasaki mode function /is(fc,?7) = C,\JlK,t\l c^ (C being the 
comoving curvature perturbation) satisfies the modified 
equation of motion 



as K 



Ms 



2,2 {a^lcs)'' 



ctk 



a^/cs 



Ms = 0. 



(5) 



Derivatives are with respect to conformal time 77, a(vi) 
is the FLRW scale factor and ei = — dlnTJ/dlna is the 
first Hubble flow function [H being the Hubble parame- 
ter) . In addition to the induced modification of the per- 
turbations' effective potential, Cg also appears as their 
propagation speed. Let us notice that for DBI inflation, 
Eqs. © and dH) give 



'P 



a?T{ip) 



(6) 



and therefore ensure that Cg < 1. In the following, we 
will only consider the k-inflation models for which Cg < 1 
(see, however, Refs. [11 [H H [H, [H, [13 ) • 

In this paper, we are interested in constraining the 
shape of the scalar and tensor primordial power spectra 
for all the kinetically modified models admitting an ob- 
servable "slowly-evolving" phase. To this end, we assume 
the existence of a background solution supporting a long 
enough inflationary era to solve the FLRW problems. 

The fact that Cg is time-dependent and appears in front 
of the k^ term in Eq. ([5]) has important consequences and 
brings new complications into the problem. In particular, 
one can no longer use the standard techniques consisting 
in solving the mode equation in terms of Bessel func- 
tions. In the case of k-inflation, one is therefore forced 
to rely on a new method of approximation. This is done 
in Ref. [l^l and in Ref. [19], this last article making use 
of the Wentzel-Kramers-Brillouin (WKB) /uniform ap- 
proximation [20, |2l|, [24I. In this manner, one can solve 
Eq. (O perturbatively in terms of a double hierarchy of 
functions 0, El [H, [M H deflned by 
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(7) 



din a 

The ti are the usual Hubble flow functions introduced 
in Refs. [2^, |2^, [27|, while the 6i are their equivalent 
based on the sound speed and therefore encode its rate 
of change. At first order in these parameters, the scalar 
and tensor primordial power spectra^ can be expressed 



^ For single field models, (^ is a conserved quantity on "super-sonic" 
length scales. As in the standard case 12811 . this allows us to 
propagate the primordial spectrum from the end of inflation to 
the time of decoupling. 
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(8) 
The wavenumber /c^ is a pivot scale around which these 
power spectra are expanded, see Ref. [19]. All the func- 
tions appearing in this equation have to be evaluated at 
the conformal time 77* at which this pivot mode crossed 
the sound horizon, namely the solution of 



Hj^ '/^ 



Csiv*) 



(9) 



Notice that although the tensor perturbations are not af- 
fected by the non- minimal kinetic terms in the scalar sec- 
tor, Cg nevertheless appears in Vh due to the evaluation 
of the flow functions at the sound horizon crossing time 
instead of Hubble crossing. In Eq. ([S]), D — 1/3 — ln3~ 
—0.76 is a constant coming from the uniform approxi- 
mation, and an overall constant factor of 18exp(— 3) has 
been absorbed into H^^f |2l|. From Eq. ([8]) we can read 
off the spectral indices to first order. 



n - 1 = -2ei - 62 + (5i , n = -2ei, 



(10) 



and the k-inflation tensor to scalar ratio r — IGeiCg. 

As expected, the power spectra in Eq. ([5]) are not the 
same than those obtained for a canonically normalised 
scalar field due to the presence of the two sound flow 
parameters Cg and Si. Notice that in the limit Cg = 1 — 
0{ei), one has 61 = C(e^) and the above spectra only 
differ from their usual analogue at second order in the 
Hubble flow functions [29]. 

In the following, wc compare the induced Cosmic Mi- 
crowave Background (CMB) anisotropies seeded by these 
power spectra to the WMAP5 data [13, [Mj, Isl IH, 
lii, 35]. Using Markov-Chains-Monte-Carlo (MCMC) 
methods, we extract the posterior marginalised proba- 
bility distributions on the sound and Hubble flow pa- 
rameters and discuss the physical consequences. More- 
over, as shown in Refs. [2J, [Sy, [37|, [S^ , kinetically modi- 
fied inflation models may induce a large amount of non- 
gaussianities in the statistics of the CMB anisotropies, 
usually measured in terms of the /^^ parameter '39'| . The 
expected value for /^^^ strongly depends on the model at 
hand: as shown in Ref. [2^, f^^ for equilateral configura- 
tions demands the knowledge of the third X-derivative of 
P{X, (f) in addition to the sound and Hubble flow func- 
tions. For these reasons, we have chosen to split our 
analysis into two parts. In the first section, we do not 
consider the non-gaussianity constraints, and our results 
are only based on the observed two-point functions of the 
CMB temperature and polarisation anisotropies. In the 



second part, we implement the WMAP5 bounds on f^^^ 
as prior knowledge on the sound speed values in the re- 
stricted class of DBI models: for those, we indeed have 
the universal expression [2J, |40, l4l| 

From now on, fj^^ is implicitly assumed to refer to 
equilateral configurations, the "local" estimator being 
/nT'^' ^ ^i^ij^i) ^^^ '^^11 below the current detection 
thresholds, see Ref. [li,!!!,!!^. 



II. CMB POWER SPECTRA 

We have used a modified version of the CAMB code J42 
to compute the CMB temperature and polarisation 
anisotropics seeded by the scalar and tensor primordial 
power spectra of Eqs. ([5]). Comparison with the WMAP 
data is then performed by exploring the full parameter 
space, including the standard cosmological parameters, 
through MCMC Bayesian methods implemented in the 
COSMOMC package [H, Hi]. The likeHhood estimator is 
provided by the WMAP team ^. 



A. Cosmological model 

The assumed cosmological model is a flat ACDM uni- 
verse involving four cosmological base parameters (with 
priors chosen from the posteriors of the previous CMB 
experiments): the number density of baryons Qh, of cold 
dark matter fldm, the reionization optical depth t, and 
9 which measures the ratio of the acoustic horizon to the 
angular diameter distance (see Ref. [43 )■ As performed 
by the WMAP team, we have now included the lensing of 
the CMB power spectra and performed a marginalisation 
over the Sunyaev-Zel'dovich amplitude Agz , in unit of the 
analytical model of Ref. 45j and with a flat prior in [0, 2] . 
For the sake of clarity, and to allow comparison with pre- 
vious works [4y|, we have, at this stage, not considered 
more data sets than the WMAP5 data, the Hubble Space 
Telescope (HST) constraint (Hq = 72±8km/s/Mpc|43|) 
and a top hat prior on the age of the universe between 
lOGyrs and 20Gyrs. Notice that we have chosen the 
pivot scale to be in the middle of the observable range, 
namely K = 0.05 Mpc"^ [H]. 



B. Observable primordial parameters and priors 

It is important to notice that the new parameters ap- 
pearing in the primordial power spectra of Eqs. ([S]) in- 
duce new degeneracies. In particular, and contrary to 
the usual case, ei is no longer uniquely linked to the ten- 
sor to scalar ratio due to the presence of Cg. Therefore, 
the product eiCg is a better motivated observable param- 
eter in k-inflation. Its order of magnitude being a priori 






FIG. 1: One dimensional marginalised posterior probability 
distribution (solid) and mean likelihoods (dotted) for the base 
and derived cosmological parameters, given the WMAP5 data 
and the k-inflation power spectra of Eq. © ■ 



unknown, we have chosen a uniform logarithmic prior^: 
— 10 < log(eiCs) < —0.5 (Jeffreys' prior). 

Apart from the overall amplitude, 'Pf(fc) in Eqs. ([5]) 
has exactly the same functional form than the canoni- 
cal scalar slow-roll power spectrum at first order upon 
the replacement ei — > ei — (5i, and £2 ^-^ £2 + <5i. In 
particular, this is the case for the spectral index Ug in 
Eq. pil)l . As a result, ei, €2 and 61 are fully degenerated 
and there are only two observable parameters, namely 
(ei ~ Si) and (£2 + Si), for which we have chosen a uni- 
form prior in [—0.3,0.3]. This important point already 
shows that, at the scalar perturbation level, there is no 
chance of constraining Si separately, neither ei nor £2; 
this is si gnif icantly different from the usual slow-roll con- 
straints [4^, I43, |50, |5l|- However, considering the pri- 



^ Recall that we are considering the case c^ < 1. 
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FIG. 3: Two dimensional marginalised posteriors (point den- 
sity) together with the one and two-sigma contours. Although 
(ei —(Si) and (e2 +5i) are not bounded, they are strongly cor- 
related to give acceptable values of the spectral index (left 
panel). The right panel exhibits the correlations between 
log(eiCs) and n^=\~ 2(ei — (5i) - (e2 + 5i). 



FIG. 2: One dimensional marginalised posterior probability 
distribution (solid) and mean likelihood (dotted) for the sam- 
pled sound and Hubble flow parameters given the WMAP5 
data. Only log(eiCs) is upper limited, this parameter fixing 
the tensor-to-scalar ratio. Notice in particular that ei is not 
constrained. The parameters (ei — (5i) and (e2 + <5i) are, how- 
ever, correlated (see Fig. [3]). 



mordial tensor power spectrum Vhik), one sees that ei 
is still uniquely linked with the tensor spectral index n,^ 
[see Eq. (fTO|) ]. As one may expect, it is unconstrained by 
the data at present, and so should be ei. We have how- 
ever included it for completeness as the third primordial 
parameter in the MCMC exploration with a logarithmic 
prior in — 5 < logei < —0.5. Finally, the overall ampli- 
tude of the CMB anisotropics is fixed by P^ — 'P(;{k^) for 
which we have chosen the usual uniform prior distribu- 
tion on ln(10i°P,) in [2.7,4]. 



C. Results 

The MCMC exploration has been stopped according 
to the R~statistics implemented in CDSMDMC [H, \E^ as 
soon as the variances between the different chains agree 
at a few percent. The chains contain 800 000 samples 
and lead to the marginalised posterior probability distri- 
butions represented in Fig. [T] to Fig. [31 The base cosmo- 
logical parameters associated with the ACDM universe 
are compatible with the recent results (S^ , although our 
primordial power spectra are different. 

Fig. [5] shows the posterior probability distributions 
(solid curves) for each of the primordial parameters we 
are interested in. As we expected, there is no bound on 
logei and the posterior is just given by the prior. On 
the other hand, the maximum amount of tensor modes 



in the CMB power spectra gives an upper limit on the 
parameter log(eiCs). Its posterior is indeed upper limited 
and at 95% of confidence, one finds 



log(eiCs)<-2.3. 



(12) 



Notice that the above upper limit is stronger than the 
ones reported in the recent literature [33, [53| . This ef- 
fect comes from our Jeffreys' prior choice on eiCg, which 
favours low values of the parameter [4J|, and possibly 
from the new natural choice of the sampling parame- 
ters (ei — Si) and (e2 + <^i)- Let us mention that this 
prior choice is motivated by our ignorance on the scale 
at which tensor perturbations contribute to the CMB, 
a priori'^. The slight discrepancies appearing in be- 
tween the posterior and the mean likelihood for this 
parameter signal some correlations with other parame- 
ters [H, Hi]- Also, both (ei - Si) and (e2 + i5i) exhibit 
unbounded distributions, and their mean likelihood (dot- 
ted curves) differs from the marginalised posterior. The 
two-dimensional marginalised probability distribution in 
the plane (ei — Si, 62 + Si) is plotted in Fig. [3l and clearly 
exhibits the above-mentioned correlations. The physical 
interpretation is straightforward: the "correlation strip" 
is given by the allowed values of the scalar spectral index 
ng [see Eq. (fTOj) ]. The barely visible thickening of this 
strip comes from some slight correlations with the pa- 
rameter P*, and then Agz, due to the functional form of 
the power spectra. Finally, as can be seen in Fig.[3l large 
values of log(eiCs) are also degenerated with the spectral 
index n^ in a way reminiscent of the corresponding corre- 
lation between tensor-to-scalar ratio and spectral index 
in standard inflation. 



^ From a flat prior on eiCs, we obtain eiCs < 0.023 at two-sigma 
level (by importance sampling). 




0.9 0,92 0.94 0.96 0.98 1 1.02 -8 -6 -4 -2 

1-2(8, -5,)- (e^ + 5,) ln[10=H/nip,] 



FIG. 4: Marginalised posterior probability distribution for 
the derived primordial parameters rig and Hint obtained by 
importance sampling. 



In Fig. m we have used importance sampling to 
reconstruct the probability distribution of n^ and 
ln(10^-ffi„f/7TT.p[). The probability distribution for n^ 
ends up being compatible with the one found by assum- 
ing a power law primordial spectrum J35i |. At 95% of 
confidence, one gets 



0.003 < 2(ei - Si) + (€2 + 6i) < 0.075. 



(13) 



This limit favours a "red tilted" scalar power spectrum, 
at two-sigma level, in presence of tensor modes and vary- 
ing sound speed*. The probability distribution associ- 
ated with iJjnf is only upper limited since this parameter 
is fully degenerated with eiCg in fixing the amplitude of 
the primordial scalar power spectrum [see Eq. ([5])]. The 
prior distribution on iJinf is logarithmically uniform and 
induced by our Jeffreys' prior on eiCg and P*. This is con- 
venient since the energy scale of inflation is also unknown 
a priori. We find 



In 10^- 



.H;, 



< -0.59, 



(14) 



at 95% of confidence. 
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FIG. 5: Marginalised posterior probability distribution for the 
DBI primordial parameters by including the non-gaussianity 
bounds on fj^^. The first Hubble flow parameter ei is now 
upper limited. 



A. New priors 

A straightforward implementation of the non- 
gaussianity bound is to include it as prior knowledge on 
the acceptable Cg values. In fact, this occurs via the so- 
called "Lorentz factor" 7^ = 1/Cg which enters into the 
expression (jlip . Since its order of magnitude is now fixed, 
we have chosen to sample the parameter space over 7^, 
instead of log(eiCs), with an uniform prior over [1,467]. 
All the other sampling parameters and priors remain the 
same as in Sect. [TTl as well as the cosmological model. 



III. INCLUDING NON-GAUSSIANITY 

In this section, we use the two-sigma WMAP5 bound 
on the /j^j^ parameter associated with "equilateral con- 
figurations" in Fourier space |35j . 



151 < /,, < 253. 



(15) 



As already mentioned, the predicted value of f^^ for k- 
inflation models is not universal and depends on the ex- 
plicit form of the third derivative of P{X, ip) with respect 
to A. As a result, we have chosen to limit our analysis 
to the sub-class of DBI models: for those, Eq. pT|) can 
be used. 



* This is, however, prior dependent; from a fiat prior on eiCs, the 
bounds become [-0.024,0.063]. 



B. Results 

The marginalised posterior probability distributions 
associated with the base and derived cosmological param- 
eters are essentially the same as in the previous section 
and we do not discuss them any further (see Fig. [T]) . The 
posteriors of the parameters appearing in the primordial 
power spectra are plotted in Fig. [5l 

In fact, the non-gaussianity prior imposed on 7^ breaks 
the degeneracies we encountered previously between ei 
and Cg, and which motivated the use of log(eiCs) as a 
sampling parameter. The sound speed being now limited 
from below, the tensor-to-scalar ratio becomes sensitive 
to the first Hubble flow parameter and we find a "weak" 
upper bound on logei. At 95% of confidence. 



logei < -1.1 



(16) 
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FIG. 6: One and two-sigma contours of the two dimensional 
marginalised posteriors for correlated pairs of primordial pa- 
rameters. The shading traces the two-dimensional mean 
likelihood. The right panel shows that the non-gaussianity 
bounds on f^^ induces an upper bound on the first Hubble 
flow parameter ei. 



On the other hand, (ei — ^i) and (e2 + Si) are still un- 
bounded but correlated along the acceptable values of 
the spectral index (see Fig. ^. As a result, Si remains 
unconstrained, and so does £2. The posteriors for the 
spectral index and the energy scale of inflation iJinf ap- 
pear to be relatively insensitive to the inclusion of the 
non-gaussianity bound. Their respective two-sigma lim- 
its are still given by Eq. ^^ and Eq. HH). 

In addition to imposing strong bounds on Cg for DBI 
models, the non-gaussianity limits of Eq. (fT5|) allow dis- 
ambiguation of the first Hubble flow parameter, and thus, 
some inference on the evolution of the Hubble parameter 
during inflation. It is clear that the limit of Eq. (|16p is 
not very tight, and quite close to the values acceptable 
for the "slowly-evolving" limit. However, improving the 
bounds on /^^^ "^iH directly improve the determination 
of the acceptable log(ei) values: one may therefore ex- 
pect tighter constraints from the next generation of CMB 
experiments. 



IV. CONCLUSION 

We have used the WMAP fifth year data to constrain 
all single field k-inflationary models admitting a phase 
of "slow-rolling" evolution during the generation of the 
cosmological perturbations. 

The probability distributions of the Hubble and sound 
flow parameters, e^ and Si, evaluated when the scalar 
perturbations crossed the sound horizon, are significantly 
different than the ones associated with canonically nor- 
malised single field models. The presence of two addi- 
tional parameters in the primordial power spectra, Cg and 
i5i, leads to new degeneracies in the parameter space: 
there is no longer any constraint on ci and £2 separately. 



neither on Cg or ^i alone. The energy scale of inflation 
is nevertheless upper limited: H■■^^i/m-p^ < 6 x 10^^ at 
95% of confidence. In fact, only two particular combi- 
nations of the primordial parameters end up being con- 
strained by the CMB power spectra: the quantity eiCg 
in Eq. (fT^ . which encodes the tensor-to-scalar ratio in 
the CMB spectra, and 2{ei — Si) + (£2 + Si), which fixes 
the scalar spectral index in Eq. (J13p . Under our prior 
assumptions, Eq. p3)) shows that the data favour a "red- 
tilted" scalar power spectrum, at two-sigma level, with 
tensor modes and varying sound speed. 

Implementing the current WMAP5 bound on non- 
gaussianity through the /^^ parameter (for equilateral 
configurations) breaks some of the above-mentioned de- 
generacies. For the DBI inflationary models, more than 
restricting the acceptable values of the sound speed Cg, 
we find that the non-gaussianity bound induces an upper 
limit on the first Hubble fiow parameter: log(£i) < —1.1 
at two-sigma confidence level. 

The present results could be improved by various 
means. A first step may be to add more data sets, but 
according to the previous discussion, this will certainly 
improve only the limits on the spectral index and tensor- 
to-scalar ratio combinations without breaking the degen- 
eracies between the ei and Si. Therefore, it is clearly 
more promising to improve the tests that bring new in- 
formation on non-gaussianities. We have implemented 
them for DBI models only through a prior knowledge 
over Cg. Clearly, our approach could be improved by di- 
rectly evaluating the three-point functions and compar- 
ing it to the observed one, provided one can define an 
tractable likelihood for MCMC exploration. Such an ap- 
proach is currently technically difficult. Evaluating the 
three-point functions from CMB data is not trivial (see 
Appendix A in Ref. [3a|); and, from a theoretical point of 
view, new parameters have to be introduced to quantify 
the non-gaussian signals in a generic way for k-inflation. 
Therefore, one faces the risk to add even more degenera- 
cies. On the other hand, with the incoming fiow of data 
on cosmological observables, this represents an interest- 
ing challenge for the forthcoming years. 
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